Search for Long-lived States in Antiprotonic Lithium 

J. Revai 

Research Institute for Particle and Nuclear Physics 
H-1525 Budapest, P.O.B. 49, Hungary 

V. B. Belyaev 
Bogolyubov Laboratory of Theoretical Physics 
Joint Institute for Nuclear Research 
I4198O Dubna, Russia 
(Dated: February 2, 2008) 

Abstract 

The spectrum of the {Li^'^ + p + 2e) four-body system was calculated in an adiabatic approach. 
The two-electron energies were approximated by a sum of two single-electron effective charge two- 
center energies as suggested in [^]. While the structure of the spectrum does not exclude the 
existence of long-lived states, their experimental observability is still to be clarified. 
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I. INTRODUCTION 



One of the most impressive success-stories of the last decade in few-body physics are 
the high-precision experimental and theoretical studies of long-lived states in antiprotonic 
helium |1|, H, || . In view of this fact it is natural to pose the question, whether such long- 
lived antiprotonic states can exist in other systems, too. There were some experimental 
attempts to observe delayed components in annihilation products after stopping of slow 
antiprotons in different media 0, however, except for helium, no evidence of long-lived 
states was found. Theoretical predictions concerning the possible existence of such states 
could largely increase the willingness of experimentalists to go on with their attempts. 

There are two possible directions in which further candidates for long-lived antiprotonic 
states could be searched for. First, one could consider atoms with electron structure similar 
to helium, that is, noble gases with closed outermost electron shells. This possibility was 
examined in some detail in |^ with an essentially negative answer concerning the possibility 
of formation of long-lived states after antiproton capture in noble gases. 

The second possibility is to consider the next simplest atom, the lithium, which has three 
electrons, one of which could be replaced by the antiproton. In order to get an idea of the 
possibility of the occurrence of long-lived states in this four-body system, we have performed 
a semi-quantitative calculation of its level structure. 

II. CALCULATION METHOD 

The Hamiltonian of the {Li^~^ + p + 2e) four-body system can be written as 




(1) 



with the two-electron Hamiltonian 



U^\r,, ra; R) = h^'\r,; R) + /i«(r2; R) + 
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(2) 



ri - r2 



The single-electron two-center Hamiltonian h^^^ corresponds to the electron motion in the 
field of Li^'^ and p fixed at a distance R: 



In eqs.(l-3) R is the vector pointing from Li^+ to p, while the rj are the electron coordinates 
measured from the Lz^+ — p center of mass. ^ M is the reduced mass of Li^+ and p: 

1 _ 1 1 

M 171^3+ rup^ 

while a is defined as 

Trip „ 
a— ^ 



mLi3+ + nip 

To calculate the spectrum of this four-body system or, at least, a part of it, we have used 
a Born-Oppenheimer-like (BO) approximation, in which the solution of the Schrodinger- 
equation 

{H-Ei)^i{v,,v-,,K)^Q (4) 
is attempted in two successive steps. First, the equation 

(ri, r^; R) - £g(i?))<l>„^(ri, v^, R) = 0, (5) 

describing the motion of two electrons in the field of Li^+ and p separated by a fixed vector 
R, has to be solved. The solutions are characterized by the conserved quantum number 

IX - the sum of the electron angular momentum projections on the R direction. Next, the 
total wave function ^;^(ri, r2, R) is approximated as a single product: 

u-^ (R) 

^;:(ri,r2,R) = ^^I)^M„-^(0,e,O)$„^(ri,r2;R), (6) 

where the Wigner's /^-functions are needed to ensure correct angular momentum quantum 
numbers JMj for the total wave function. The energy eigenvalues E^^ are calculated from 
the radial equation 



+ <(i?) = 0, (7) 



which is obtained by substituting Eq.(6) into Eq.(4), multiplying by 
G, 0)*$„^(ri, R)* and integrating over the electron coordinates and the 
angular variables {(f), ©) of R. In deriving Eq.(7), according to the adiabatic approximation. 



^ These are not exactly the Jacobian coordinates of the system, however, the extra terms in the Hamiltonian 
arising from the difference are very small and the accuracy of our calculation does not necessitate their 
consideration 

^ All masses are expressed in units of electron mass 
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terms containing the derivatives of $„At(i"i7 R) with respect to R are neglected. In BO 
type calculations usually the lowermost electron configurations are used, for which /i = 
(so called a-term) in which case the D-function in Eq.(6) reduces to a spherical harmonics 

yjMj- 

For a one-electron problem Eq.(5) is replaced by 

(/i«(r,;R) K,(r,;R)=0 (8) 

and it can be solved exactly by separation of variables in spheroidal coordinates <^n4i{T^i] R) 
being the well-known two-center wave functions. In our case, however, the problem of two 
interacting electrons in the field of two fixed Coulomb-centers is not solvable exactly and the 
determination of eigenvalues and eigenfunctions necessitates a highly non-trivial calculation. 
One could think of a variational approach of the type [0; another possibility could be the 
diagonalization of the electron-electron interaction on the basis of two-center functions. In 
both cases the calculations are quite cumbersome and slowly converging when the size of the 
basis is increased. Since we believe, that the question of possible existence of metastability 
in antiprotonic lithium can be studied by a semi-qualitative exploration of its spectrum, we 
have chosen a simpler, elegant and efficient method of calculation for e^^j^ (R) proposed in . 

The main idea of the method is to represent the two-electron energy eigenvalue as a sum 
of two single-electron energies: 

e^'\R)=e^zUiR) + e^^UR)^' (9) 

where e^l Z2 (R) ^^z^z2 (R) energy eigenvalues of two-center equations of the type 
(8) with effective charges (^1,^2) and (2:1,-22) instead of the physical charges (3,-1). The 
effective charges (^1,-22) are chosen in such a way, that €^^^^^{R = 0) and e^^^^^{R = 00) 
should reproduce the experimental values of the first ionization potentials of He atom and 
Li~^ ion, respectively. Thus we get zi = 2.3578 and Z2 = —1.0135. As for (^1,^2), the 
corresponding two-center eigenvalues in the R = and R = 00 limits should reproduce the 
second ionization potentials of He and Li~^, and the second electron is in this case the last 
one, therefore the physical values Zi = 3 and Z2 = — 1 were taken. The electron-electron 
repulsion is taken into account in this method by the deviation of the effective charges {zi, Z2) 

^ Here and in the following, where no confusion can arise, the indices {n^) arc omitted 
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FIG. 1: Electronic energies of the {Lr"^ + p + 2e) and (Li^+ + p + e) systems 



from their physical (integer) values. In this way the approximate e^^-* [K] of Eq. (9) reproduces 
the experimental two-electron binding for the two limiting cases = and R = oo, while 
for intermediate i?-values the solution of the corresponding (effective) two-center problems 
seems to provide a reasonable interpolation prescription. This approach has been checked 
in the case of {He^^ + p + 2e) and {Lr"^ + p + 2e) systems; its results were compared with 
those of a detailed variational calculation of [0. Results of the calculations are presented 
in Fig.(l). The agreement of e^'^\R) obtained from Eq. (9) with the variational values is 
amazingly good in a wide range of R. The same procedure was applied to calculate the 
energy of the first excited electron configuration, where the limiting cases were adjusted to 
reproduce the energies of the first excited (ls2s) states of He and Lz+. 

Having obtained the electronic energies e'^'^\R) the effective potentials of Eq.(7) 



Jn^ _ J{J+^) _ 1 , ^(2) / ^>^ 



(10) 



can be calculated. Here the index /i = was omitted, while the electron configuration label 
n can take the values n = (Isls), (ls2s). For some values of the total angular momentum 
J the effective potentials are shown in Fig. (2) for the ground- and first excited electron 
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FIG. 2: Effective potentials in the (Li'^^ + p + 2e) system 



configuration. The energy eigenvalues E;^" are then calculated by solving Eq.(7) with these 
effective potentials, and the "vibrational" quantum number v is introduced to distinguish 
among the states with the same J-value. 

III. RESULTS AND DISCUSSION 

The resulting spectrum of the {Li^~^ +p + 2e) system is shown in Fig. (3). Apart from 
the energy levels of the initial system (full black circles) we have shown also the energies of 
daughter states which can be formed after Auger-emission of one or two electrons. The open 
circles correspond to the energies of the (Li^^ +p + e) three-particle system, while the open 
squares are the hydrogen-like two-particle energies of the {Li^~^ +p) system. Comparing the 
spectra of Fig. (3) with the well-known spectra of the He atomcules (see e.g. 0), we can 
find some apparent similarities and differences. The basic similarity can be formulated as 
follows: there are many states in the spectrum, from which Auger-emission is possible only 
with large electron orbital momentum and therefore is strongly suppressed. This could be 
one reason for metastability of these states; of course, this is only a necessary condition and 
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FIG. 3: Spectra of antiprotonic lithium atoms 
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by no means a sufficient one. 

Tfie basic difference, on tfie otfier fiand, is tfie mudi fiiglier density of states in tfie 
expected capture region (around tlie Li atom ground state energy) wfiidi is due to tfie 
essentiai difference in the electron structure of He and Li: the last efectron is strongiy 
bound in He, while very loosely in Li. 

It can be noted, that the spectrum of the {Li^~^ + p + e) system (open circles) strongly 
resembles the He atomcule spectrum, therefore long-lived states in an isolated {Li^'^ +p + e) 
system could be certainly expected. However, in contrast to the He case, this system is 
charged and thus its interaction with atoms of the surrounding medium might be more 
violent, leading to a faster coUisional de-excitation of these states. 

IV. CONCLUSIONS 

We have calculated the spectra of (Li^+ +p + 2e) and (Li^"*" +p + e) four- and three-body 
systems. Although the structure of the obtained spectra allows the existence of long-hved 
states, our calculations do not put us in a position to make definite statements about their 
experimental observabihty. This latter depends on several further factors, as well. One is 

the formation mechanism: for the time being we have no information about the population 
rate of the huge amount of states in the vicinity of expected antiproton capture energy. The 
physics of formation of the {Li^~^ + p + 2e) system in the reaction 

(Li^+ + 3e) +p — > {Li^^ +p+2e) + e 

is quite different from the analogous process in He due to the large difference in the binding 
energies of the outermost electron. In the Li atom the first ionization potential is only 
0.198 a.u., which means that even adiabatic ionization is possible: when the distance between 
the antiproton and Li atom becomes less than 5 a.u., the binding energy of the last electron 
in their common field becomes zero and the electron is "pushed" into a continuum state. 

Another unknown factor is the way, how the eventually formed (Li^+ + p + 2e) systems 
interact with the media atoms and what is the role of collisional de-excitation in their 
life-time. In any case, in order to reduce the undesired effect of this factor, probably, the 
experiments looking for long-lived states should be performed in dilute vapors of Li. 
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